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Abstract 

We study the anapole moment of the lightest neutralino in the 
constrained Minimal Supersymmetric Standard Model (cMSSM). The 
electromagnetic anapole is the only allowed electromagnetic form fac¬ 
tor for Majorana fermions, such as the neutralino. Since the neutralino 
is the ESP in many versions of the MSSM and therefore a candidate 
for dark matter, its characterization through its electromagnetic prop¬ 
erties is important both for particle physics and for cosmology. We 
perform a scan in the parameter space of the cMSSM and hnd that 
the anapole moment is different from zero albeit very small (< 10“^ 
GeV“^). Combined with experimental constraints like the Higgs mass 
and the DM relic density, the allowed region of parameter space lies 
within the reach of future direct DM searches. Thus, the anapole mo¬ 
ment could be used as a complementary constraint when studying the 
parameter space of the cMSSM and other similar models. 


1 Introduction 

One of the best motivated extensions of the Standard Model (SM) is the 
Minimal Supersymmetric Standard Model (MSSM), since, besides giving a 
solution to the hierarchy problem, it provides us with a good candidate for 
cold dark matter (CDM), namely, the lightest neutralino. 
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A number of different experiments are working (or will be soon) in the 
search for a direct or an indirect signal of dark matter (DM) (for recent 
reviews on dark matter detection see 0|)' If DM is detected, there will 
be a need to differentiate between the candidates, characterizing them as 
much as possible. In the last few years, there has been intense work on the 
electroweak properties of dark matter since they might be relevant for the 
calculation of DM decays and annihilations, (3 15 , which have consequences 
in astrophysical processes 16,17 and are important in indirect astrophysical 
searches for DM, as in the calculation of the annihilation cross section of the 
DM itself. 

Motivated by this, since 2009 we have been studying the toroidal dipole 
moment (TDM) of Majorana particles, which is related to the anapole mo¬ 
ment, one of the least studied electromagnetic properties of a particle 


18-21 


Lately, there has been a surge of interest in the study of anapole moments 
from the astrophysical as well as the particle physics points of view (e.g. 


22 |^). 

Since hrst discussed in 1957, the anapole moment, introduced by Zel’dovich 
has been investigated in different helds of science and technology. The 


26 


anapole moment corresponds to a T invariant interaction, which is C and 
P non-invariant [^. Within the Standard Model it has been calculated in 
neutrino and hadron physics 27 30 . In nuclear physics it has been studied 


in atomic nuclei 31 33 . In relation to the DM problem, there are different 
proposals to observe this “physical observable”, extracting its value from di¬ 
rect measurements between DM and atomic nuclei 23,24,34,35 . Finally, in 


engineering various aplications have been studied, in areas such as “Ferrite 
electromagnetic properties in dielectric nanoparticles 


resonators” 36 
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and in the study of electromagnetic radiation in antennae with a helicoidal 
toroidal geometric distribution 


38 -42 , among others. 


Recently, Ho and Scherrer have proposed that dark matter interacts with 
ordinary matter exclusively through the anapole moment j^. They calcu¬ 
late the anapole moment needed to obtain the right amount of DM relic 
abundance, and the anapole DM signatures that could be observed in the 
LHC 43 . Haish and Kahlhoefer have shown the importance of loop contri¬ 


butions to the scattering cross section of dark matter, in particular those in¬ 
duced by the anapole interaction 44 . More recently, del Nobile et al made a 


halo-independent analysis of direct DM detection data considering that it has 
only anapole and magnetic moment dipole interactions 24 . Also, an analysis 


on the loop corrections for leptophilic DM and internal bremsstrahlung was 
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presented in 23 , where the authors calculate the DM anapole and dipole 


moments in a toy model using direct detection data. 

In this paper we calculate the anapole moment of the neutralino at the 
one-loop level within the constrained Minimal Supersymmetric Standard 
Model (cMSSM). In the MSSM the anapole contributions of the neutralino 
arise exclusively through radiative corrections to the vertex XXI- We do a 
scan in the five parameter space of the cMSSM and compare the results for 
the anapole moment with the above mentioned experimental limit. In our 
analysis we take into account also other experimental constraints, namely the 
Higgs boson mass and the decays 6 —)■ sy and B —)■ in order to find the 

viable regions of parameter space. We hnd that although the anapole mo¬ 
ment is very small througout the regions studied, it is possible to distinguish 
between the different regions of parameter space through it, which makes it 
indeed an important property when characterising dark matter. Our results 
agree qualitatively with those found by Ho and Scherrer 


22 


The article is organised as follows: in section H we present a very brief 
summary of some aspects of the constrained MSSM relevant to our calcu¬ 
lation. In section HI we review the general form for the electromagnetic 
vertex of a particle, and in particular for a Majorana particle. We introduce 
the anapole moment and its relation to the toroidal dipole moment. In sec¬ 
tion IV we explain the methodology used to calculate the anapole moment 
of the neutralino in the cMSSM and evaluate it for different values of the 
parameters. Section V presents the obtained results and our conclusions. 

2 The MSSM and the neutralino as candi¬ 
date for dark matter 


According to the latest results of WMAP, the CDM density is 45 

noMh^ ~ 0.1109 , (1) 

where h is the Hubble constant in units of 100 km sec“^ Mpc“^. The ther¬ 
mally averaged effective cross section times the relative speed of the dark 


matter particle, needed to get this relic density is 46-48 

< av >oc 


( 2 ) 


consistent with the assumption of a weakly interacting dark matter particle 
(WIMP) with a mass between 10 GeV - (few) TeV. 
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The minimal supersymmetric extension of the Standard Model (MSSM) 
provides us with one of the best WIMP candidates for dark matter: the ligh- 
est neutralino (for reviews on SUSY see for instance 49,^). The MSSM 


requires two complex Higgs electroweak doublets to give mass to the up and 
down type quarks in order to avoid chiral anomalies. The MSSM has also 


a new discrete symmetry, R parity, dehned as R = (—1) 


3B+2S+L 


where B 


and L are the baryonic and leptonic numbers respectively. This symme¬ 
try assignes a charge -|-1 to the SM particles and -1 to the supersymmetric 
partners, making the lightest supersymmetric particle (LSP) stable. 

However, supersymmetry has to be broken or it would have already been 
observed. To break it explicitly without the reappearance of quadratic diver¬ 
gencies, a set of super-renormalizable terms are added to the Lagrangian, the 
so-called soft breaking terms. The Lagrangian for the soft breaking terms is 
given by 

Csoft = + c.c. - , (3) 


where Ma are the gaugino masses, and are trilinear and bilinear cou¬ 
plings respectively, and (m^)® are scalar squared-mass terms. It is assumed 
that supersymmetry breaking happens in a hidden sector, which communi¬ 
cates to the observable one only through gravitational interactions, and that 
the gauge interactions unihy. This means that at the GUT scale the soft 
breaking terms are “universal”, i.e., the gauginos Ma have a common mass, 
as well as the scalars (m^)® and the trilinear couplings, Requiring elec¬ 
troweak symmetry breaking hxes the value of R®-^ and the absolute value of 
the Higgsino mixing parameter \fi\. This is known as the constrained MSSM 
(cMSSM) which is described by hve parameters: the unihed gaugino mass 
mi/ 2 , the universal scalar mass mo, the value of the universal trilinear cou¬ 
pling Ho, the sign of Higgsino mass parameter /i, and the ratio of the vacuum 
expectation values of the two Higgses, tan jS. 

After the electroweak symmetry breaking the neutral and charged states 
in the MSSM can mix. In the case of the neutral ones they give rise to 
a set of four mass eigenstates, the neutralinos. It is the lightest one of 
these that is the LSP and a good candidate to dark matter in many SUSY 
models. The lightest neutralino, in the gauge eigenstate basis, is a function 
of the neutral higgsinos and the neutral gauginos (wino and bino) and its 
properties will depend on the mixing, which in turn depends on the soft 
breaking parameters. 
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Before WMAP, the cMSSM was compatible with the limit ~ 

0.1 — 0.3 and other direct and indirect low energy and collider data in a hnge 
parameter space region called the “bnlk”. However, after the constraint by 
WMAP to VL£)m fib'll and with the recent limits to the sparticles masses from 
the LHC, which excludes light masses, the bulk region in tuq — mi /2 is no 
longer viable. Moreover, in the cMSSM the LSP neutralino is, practically for 
all cases, an almost pure bino state which annihilates itself more efficiently 
into leptons through right hand sleptons due to their higher hypercharge. 
However, with the newest data from WMAP, this mechanism is not suffi¬ 
ciently efficient. There are still three favoured scenarios that require some 
very specific accidental relations between some parameters at the electroweak 
scale. 

In the cMSSM at low mo, there is a region with almost degenerate r —Xi- 
In this case the populations of these two particles are almost the same, mak¬ 
ing the NLSP f thermally accesible. The mass diference between the scalar 
tan and the lighest neutralino, AM = m^—m^o, controls the population ratio 
of these two species through the Boltzmann factor exp(—AM/T/). There¬ 
fore, it is a very sensitive parameter which enters into the calculation of the 
relic density. Whenever the coannihilation takes place, through the partici¬ 
pation of f in processes like fiXi —t ry or even fifi —> rf, the relic density 
can be reduced in comparison to the case of the bulk scenario. In this region, 
the LSP neutralino is mainly bino with a mass esentially stablished by Mi up 
to corrections of order M|//i (/x is high). The approximate formulae for the 
mass of the neutralino and the mass of fi suggest that degeneration happens 
for mo 0.145 mi/2. 

A sudden increase in the usual mechanism of coannihilation to reduce the 
relic density can occur if m^o is close to a pole. Faster and more efficient 
annihilation can take place through Higgs resonance. Given the Majorana 
nature of the neutralino, the resonant enhancement is obtained only via 
the pseudo-scalar Higgs boson. The colliders constrains to the LSP in the 
cMSSM allow the heavy “Higgs funnel”, where XiXi —t A —)■ bb/rf, which 
happens for high tan 13. Therefore, the quantities that establish this scenario 
are the quantity 2m^o — m^ and the amplitude of the pseudo-scalar, since 
they dehne the resonance prohle of A. 

In most of the cMSSM, n is very high. However, you can exceptionally 
have that fi ~ Mi, which allows much more efficient coannihilation through 
reactions as XiXi —t WW/ZZ/Zh/tt. This happens in the so-called focus 
point region where mo is very high. The focus point region corresponds 
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to high values of mg close to the border of viable electroweak symmetry 
breaking, where the value of /i decreases rapidly. When /r ~ Mi, M 2 , the 
LSP has a signihcant fraction of higgsino, and the next lightest sparticles (y® 
o Xi) have also a signihcant component of higgsino and are not much heavier 
than the LSP. Thus, the coannihilation channels are favoured. However, 
coannihilation cannot be very efficient, otherwise the relic density would 
be less than what is actually measured. In this scenario, all the sfermions 
are very heavy (more than 4 TeV) to be accesible to one of the proposed 
colliders. The LSP mass goes from close to 150 to 350 GeV, with higgsino- 
type neutralinos being 100 to 50 GeV heavier. From the perspective of a 
linear collider, an energy of more than 800 GeV is needed to reveal some of 
the properties of this scenario. The pseudo-scalar has a mass higher than 1 
TeV and very likely would not be found directly in the LHG. 

There have been several scans of the cMSSM and other MSSM type mod¬ 
els using different criteria to check its viability. Many of these are either 
Bayesian (e.g. 51-^) or frequentist (e.g. 56-^), or both ( [^|^), but 


there are also many studies using other techniques or constraints to perform 
the analysis (e.g. 61 -^). In the cMSSM it is assumed that DM is composed 
by only one type of particle, the LSP, which is usually assumed to be the neu- 
tralino. Most of the studies then look for regions of likelihood with different 
boundary conditions at the GUT scale, i.e. values for mg,mi/ 2 ,tan/3,/r and 
Aq, and phenomenological constraints at low energies, which usually include 
a combination of the following: the branching ratios b —>■ sy, B —>■ the 

dark matter relic density, the requirement of radiative electroweak symmetry 
breaking, the Higgs mass, and a solution to the g — 2 problem (although not 
all of these constraints may be considered together in every analysis), and 
constraints coming from direct and indirect searches for dark matter. From 
these various scans it is clear that the cMSSM is highly challenged, although 
there are still regions of parameter space allowed, depending on which low 
energy constraints are used. In general, to satisfy most of the above con¬ 
straints, a heavy supersymmetric spectrum is expected, with Aq ^ 0 and 
large tan (3 (for a recent overview on constraints on the cMSSM and other 
SUSY models see for instance 


56,58,66, CT) 
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3 Anapole Moment 


For 1/2-spin particles the most general expression for the electromagnetic 
vertex function, which characterizes the interaction between the particle and 
the electromagnetic held, is: 


rM(?) = /q(9^)7m + 

+ fA{,q^){,q^l^, - g'g/.)75, 


( 4 ) 


where fq^q^), f^i^q^)-, fsiq^) and fA{.q^) are the so called charge, magnetic 
dipole, electric dipole and anapole form factors, respectively; in here = 
p/ —is the transferred 4-momentum and = (i/2) [7/i, 1u] 29,68 . These 
form factors are physical observables when —)■ 0, and their combinations 
dehne the well known electric charge (Q), magnetic dipole (p), electric dipole 
(d) and anapole (a) moments. 

However, the electromagnetic properties of the neutralino (which is a 
Majorana particle) are described by a unique form factor, the anapole, fA^q^)- 
This is a consequence of CPT-invariance and the C, P, T properties of 
and the interaction Hamiltonian. Thus, the electromagnetic vertex function 
of the neutralino can be written just as 


^/.{q ) = fA{q )(g 7 m - ?'?m)75- 


( 5 ) 


The anapole moment was introduced by Zel’dovich to describe a T- 
invariant interaction that does not conserve P and C parity [^. In contrast 
to the electric and magnetic dipole moments, the anapole moment inter¬ 
acts only with external electromagnetic currents In the non- 

relativistic limit, the interaction energy with an external electromagnetic held 
takes the form 


'Hint oc —p {a ■ B) — d {a ■ E) — a {a ■ \/ X B), 


( 6 ) 


where B and E are the strength of the magnetic and electric helds, and a 
are the Pauli spin matrices. 

The anapole moment does not have a simple classical analogue, since 
/^(g^) does not correspond to a multipolar distribution. A more convenient 
quantity to describe this interaction was proposed by V. M. Dubovik and 
A. A. Cheshkov 27 : the toroidal dipole moment (TDM), T(g^). For a 


comprenhensive review on complete electromagnetic multipole expansions. 
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Figure 1: Current configuration with a toroidal dipole moment. The arrows 
on the torus indicate the direction of the current, and the TDM is directed 
towards the axis of symmetry of the torus. 


including toroidal ones, see 69 . The TDM and the anapole moment coincide 


in the case of m* = m/, i.e. the incoming and outgoing particle are the same. 
This type of static multipole moments does not produce any external helds 
in vacuum but generate a free-held (gauge invariant) potential 29 , which is 
responsible for topological effects like the Aharonov-Bohm one. 

The simplest TDM model (anapole) was given by Zel’dovich as a con¬ 
ventional solenoid rolled up in a torus and with only one poloidal current, 
see Fig. [T} For such stationary solenoid, without azimuthal components for 
the current or the electric held, there is only one magnetic azimuthal held 
diherent from zero inside the torus. 

As mentioned in the introduction, the anapole moment is a very useful 
quantity in nuclear physics, where it has been widely studied, as well as in 
astrophysics and engineering. In particle physics it is important in DM de¬ 
tection, since the DM candidates can have couplings to nuclear spins. There 
are also limits to detection of anapole dark matter in the LHC 43 , which 
exclude it for masses < 100 GeV. 

To measure the anapole moment of DM, direct detection is needed, where 
the resulting cross section from the scattering of the DM particle with a nu¬ 
cleus is measured, and from there the anapole moment or bounds to its value 
can be extracted. The hrst such upper bound was calculated in ref. 70 , using 


data from the DAMA/LIBRA collaboration 71 and from the Ge detector 













72 


of the CDMS collaboration 


4 X 10 ^fm 


(7) 


for a WIMP mass of 100 GeV (3 for the Ge detector and 4 for the Nal one). 
There are currently several experiments exploring direct DM detection 


). The best ex- 
and LUX l75 


(for a recent review on indirect and direct DM searches see 
elusion bounds for DM at present come from XENONIOO [74 
Although the anapole moment is a very small quantity, it is expected that 
the improvement in the sensitivity of future direct DM detection experiments 
will allow to put more stringent bounds on its value. In this respect, knowing 
precisely the neutron and proton spin contents of relevant nuclei is important 
for the correct interpretation of the data (see for instance 33 -35 , and refer¬ 


ences therein). It is expected that both XENONIT 76 and LUX-ZEPLIN 
(LZ) [^ will improve by ~ 100 times their measurement of the WIMP- 
nucleon cross section. In the case of XENONIT the sensitivity is expected 
to be 2 X lO”"^^ cm^ at a WIMP mass of 40 ~ 50 GeV [T^IT^, whereas LZ 


has a projected sensitivity of 10 cm^ for its full 1000 day exposure 77(79 


4 One-loop calculation 

A neutralino is a Majorana particle, and it is necessarily electrically neutral. 
This fact does not allow for a tree level electromagnetic coupling. Therefore 
the electromagnetic properties of the Majorana particle -the anapole- arise 
only via loop contributions. The anapole moment of the neutralino may 
be dehned in the one-loop approximation in the MSSM by the Feynman 
diagrams shown in Figs. and where / represents the charged fermions 
of the SM. Taking each fermionic family separately we obtain 96 Feynman 
diagrams in total, corresponding to self-energies and vertex corrections. 

We use FeynCalc [^ to calculate the amplitude of these diagrams. Since 
we are only interested in the terms that contribute to the anapole form 
factor, we isolate the ones that have the Lorentz structure 7 /^ 75 . One of the 
hrst results we obtain is that the self-energies 777 °, 7 / 1 °, 7 ^° and 7 ^° do 
not contribute to the calculation at all. If we call Sj the coefficient that 
multiplies for the ith diagram, then we have that 

= ( 8 ) 
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Figure 2 : One-loop vertex corrections to the process 7 


X°iXv 


To obtain the anapole moment a = /a(0) we use the I’Hopital rule and 


get 


a = /a( 0 ) = lim 

q ‘^^0 


E. 




q2 Qq2 


(9) 


Two- and three-point Passarino-Veltman scalar functions arise in the calcu¬ 
lation of each diagram. The two-point PV scalar function is defined as 

r d^k 


R 2 2 \ _ (27rp.)^ 

iJo[q , ^2) — 


ITT^ 


— ml] [{k -I- g )2 — ml] 


21 ’ 


( 10 ) 
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Figure 3: One-loop corrections to the self-energy for the process 7 
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and the three point PV scalar fnnction is dehned as 




(27rp) 


4-D 


d^k 


(11) 




— rn\] [{k + pi)^ — ml] [{k + ^ 2 )^ — m^ 


where A = p\,B = (pi — P 2 Y ^ — P 2 - The self-energies contain 
two point Passarino-Veit man scalar fnnctions of the type Bq {q‘^,x‘^,x^) and 
Bq (0,x^,x^). Likewise, the contribntions to the vertex corrections have two 

and three point scalar fnnctions of the type Bq , x^, x"^), Bq 


2 2 2 
mtoyy ,X^ 


2 2 2 
" , mto,m-o,x 

xr xr 


^ x‘^,y‘^ 


In both cases x and y represent the masses 


and Cq 

of the particles in the loop. 

When evalnating ([^, derivatives of the Passarino-Veltman fnnctions ap¬ 
pear. To evalnate the Bqs, as well as their derivatives, we nse LoopTools 81 


To evalnate the Cqs and their derivatives we expand them in a power series 
aronnd q^ = 0. In this way it is possible to find an analytic approxima¬ 
tion which coincides with the fnll expression in the limit = 0, simplifying 
enormously the calculation (see appendix). 

In all regions of parameter space (except for mi /2 3> mg, which is ruled 
out by cosmological constrains since the LSP is charged) it was observed 
that the four triangle diagrams involving r in the loop are almost completely 
dominant. (Fig. shows two of these diagrams. The other two are equal, but 
with the flow arrows going counterclockwise due to the Majorana nature of 
the neutralino.) The approximate analytical expressions for the contributions 
of these diagrams are 


-k 


- 4m|o 

Xi 


2m? -F 2m|o)i?o(g^, ml, ml) 2(3m?o - m?)i?o('m?o, m^, m?) (12) 


2 

x'i' 


+2 


Xi’ 

2 


(m-o-m-) -qm-o Co{q ,m^o,m-o,m^,m^,m-) - {q -4m-o) 


and 


4m?o 


xi 


2m? - 2m|o)5o(g^m?,m?) -F 2(m?o - m?)So(m?o, m^, m?) (13) 


-2(m-o-m^) Co{q , m-o, m-o,m^, m^, m^) + (g - 4m-o) [> , 
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Figure 4: Two of the four dominant Feynman diagrams for the calculation 
of the anapole moment of the neutralino. The other two are equal but with 
internal arrows directed counterclockwise. 


where k is given by 


[Zi^i (Z 2 J + Zi^ 2 ) cos 0w sin + ^ 2 , 1 ^ 1,2 cos^ 0], 

_r_ 

1 OS.7r2 r*/^c2 /y2 


O r72 


These expressions depend on the masses of the particles involved; Zjj, the 
elements of the neutralino mixing matrix; Of, the r mixing angle; as well as 
the electroweak angle 6w 

We evaluate the anapole moment within the cMSSM using Suspect 


82 


by fixing the value of Aq, tan/? and sign/r, and scanning over the other two 
parameters, mo and mi/ 2 . We then vary Aq, tan/?, and sign/r, and repeat 
the procedure. These values we then input into our own code to calculate 
the anapole moment. This code includes all diagrams contributing to the 
anapole moment, evaluated using the approximation given in the appendix 
for the Passarino-Veltman Cq functions in the limit —)■ 0. The expressions 
for each diagram are not shown explicitly in this paper. 

We have not considered in our analysis the region where m -q and rrif 

Xi 


are degenerate or quasi-degenerate, since this corresponds to a pole in the 
anapole moment function. Other methods than the one we used here should 
be employed to analyse this region. Despite the fact that sign/x > 0 may solve 
the problem of the discrepancy between the measured value of g — 2 of the 
muon and the one predicted by the SM, this does not mean negative sign/i 
is ruled out since this problem might be solved through other mechanisms, 
therefore sign/x < 0 should be also taken into consideration. 
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Figure 5: Dependence of the anapole moment (per nuclear magneton) with 
the rrir/rn^o (red) and mo/mij 2 (blue) ratios, for tan/3 = 50. We show only 
the values where the LSP is the lightest neutralino, and not the stau. The 
dependence for tan /3 = 10 is similar. 


Although the full expression for the anapole moment is a complicated 
function of the various parameters, it depends mainly on the relative values 
of and ruf, which in turn depend on the ratio of rnQ/mii 2 , as can be seen 
in Fig. 1^ No dependence on Aq or sign/i was found, however, the anapole 
moment does depend slightly on tan (3 as can be seen in Fig. 

Cosmological and experimental constrains have highly reduced the al¬ 
lowed regions of parameter space of the cMSSM. The most recent of these 
constrains is the one from the CMS and ATLAS collaborations on the mass 
of the lightest Higgs boson rrih = 125.8 ±0.6 GeV 83 -86 . This constraint is 
derived from a combination of 5.1 fb“^A/s = 7 TeV data and 12.2 fb“^\/s = 8 
TeV data. In our calculation we take this constraint as rrih = 126 ± 3, where 
the uncertainty comes from a combination of the experimental and theoreti¬ 
cal determinations of the Higgs mass. 

There is also a new measurement of 


BR{B, /x+/i") = (3.2 ± 1.5) X 10 


-9 


( 16 ) 


from the LHCb collaboration, derived from Ifb ^ of data at ^/s = 7 TeV 
collision energy and l.lfb”^ of data at ^/s = 8 TeV collision energy 87 


The excluded region due to this constraint in the cMSSM has already been 
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Figure 6: Anapole moment as a function of tuq and mi /2 (Aq = 0 and 
sign/x > 0 for tan/3 = 50 (blue) and tan/3 = 10 (green). The left panel shows 
the mo — a plane projection while the right panel shows the mi /2 — a one. 
The solid line gives the maximum value of a for that particular tan/3. 


determined in [^. We also impose the constraint coming from the branching 
ratio of 6 —)■ s'j, whose value is given by the Heavy Flavour Averaging Group 
(HFAG) as |8^ 

BR(6 ^ sy) = (3.55 ± 0.24;°;?^ ± 0.03) x 10"^ (16) 


This result can be determined directly form Suspect. Moreover, if we assume 
that neutralinos make up all of the dark matter in the universe, the WMAP 7- 
year dark matter relic abundance value = 0.1109±0.0056 45 puts even 


more strict constraints. We calculated this value using micrOMEGAs 90 


In general, the “surviving” regions have a very small value for the anapole 
moment of the lightest neutralino (10“®-10“^ GeV“^), which is consistent 
with the results of Ho and Scherrer . 

In the upper plot of Fig. we show the anapole moment values for 
different regions of parameter space in (mo,mi/ 2 ) planes in the cMSSM for 
tan/3 = 10. On top and closely around the pink dot-dashed line the stau 
and neutralino masses are degenerate, and we do not calculate the anapole 
moment in this region. The plotting program extrapolates between the values 
on both sides of the line, but there is actually a gap in the data there. 

The different phenomenological constraints are shown as follows: the re¬ 
gion above the red dashed line is where the Higgs mass is rrih = 126 ± 3 GeV, 
to the left of the pink dot-dashed line the LSP is charged, under the dotted 
blue line is the region excluded by the value of 6 —)■ sy, whereas the region 
under the dotted green line is excluded because it does not comply with the 
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Figure 7: Value of the anapole moment in (mo, ■ 07 . 1 / 2 ) planes in the cMSSM 
for tan/3 = 10 (upper plot) and tan/3 = 50 (lower plot), assuming 4o = 0 
and /i > 0, see text. The region that complies with all the phenomenological 
constraints would be to the extreme right of the plots, between the white 

lines. 
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requirement of radiative electroweak symmetry breaking. The region where 
the relic LSP density falls within the range allowed by WMAP is marked 
with a white line, while a more loose constraint, < 0.12, assuming the 
LSP is not the only component of CDM, is delimited by a white line. The 
lower plot shows the same regions, but with tan (3 = 50. 

In the two graphs we can see that the anapole moment is < (P(10“^)GeV“^ 
for every region of the space of parameters. That is the case for the bulk (low 
mo and low mi/ 2 ), already excluded by other constraints. This is also true 
for the coannihilation region (low mo and higher mi/ 2 ), where the masses of 
the and the f are almost degenerate. There seems to be a mechanism, 
which is a function of m|o — suppressing the contributions from this 
region. As can be seen from the approximate analytical formulae, there is a 
dependence on the mass difference between the stau and the neutralino both 
in the numerator and the denominator. It has to be remembered though that 
the anapole moment depends on the derivatives of these expressions, which 
obscures the mechanism at hand. The anapole moment gets relatively large 
in the focus point region, which corresponds to high mg and low mi/ 2 , close 
to the border of viable EWSB (green dotted line). The same behavior is seen 
for different tan /?, although the regions might differ in position and size. 

In Fig.j^we show the anapole moment dependence on the neutralino mass, 
and the exclusion curves of XENON and LUX 
of XENON IT andLZ 
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23 , plus the projected reach 


The regions above the lines of LUX and XENON 
are the values already excluded by the non-observation of dark matter. The 
orange points are all the ones where the lightest neutralino is the LSP, the 
brown ones are the ones that comply with the loose relic density constraint, 
and the green ones comply with the range of values we took for the Higgs 
mass. As can be seen from the graphs, the region in parameter space that 
complies with the Higgs mass constraint is at the border of the reach of LUX 
and XENON, but within the reach of LZ and XENONIT. For tan/3 = 10 
there is practically no region where there is an overlap between the loose relic 
abundance and Higgs mass constraints, so this region is basically excluded 
by all three constraints (relic abundance, Higgs mass, and lack of observation 
of dark matter). For tan/3 = 50 there is a small overlap region between the 
Higgs mass and relic density constraints. This region lies at the border of 
the exclusion region of LUX and XENON, but within the projected reach 
of LZ and XENONIT. If in the future it is possible to measure the anapole 
moment of a DM candidate, this would give us extra information on the 
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Figure 8: The plots show the dependence on the anapole moment (per nuclear 
magneton) on the neutralino mass for tan/3 = 10 and 50. The regions above 
the solid lines correspond to the excluded regions by XENON and LUX, 
and the ones above the dashed lines correspond to the projected reach of 
XENONIT and LZ. (See text for an explanation of the different regions.) 
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allowed parameter space of the cMSSM. 

The cMSSM may be too constrained to be realistic, however, using it as 
a test model, we can see that the anapole moment is indeed different for 
different regions of parameter space, and is within the reach of the future 
experiments. Thus, anapole analysis can be used as another criteria to study 
the parameter space. Although the anapole moment is insensitive to Aq 
and signp, other observables, like the Higgs mass and some decays, are not 
[^[^|^[^|^[^. This will give different exclusion regions for different 
values of the cMSSM parameters. 


5 Conclusions 


We calculated the anapole moment of the lightest neutralino in the frame¬ 
work of the cMSSM. Even though this is perhaps an unrealistically con¬ 
strained model, it is one of the most studied SUSY models, and therefore it 
is important to pass it through all possible tests. 

We found that the anapole moment of the neutralino is sensitive to mo, 
mi /2 and tan/?, but non-dependent on Aq and sign/i. The parameter space 
we scanned gives rise to an anapole moment consistent with the upper limit 


obtained by Pospelov and ter Veldhuis 70 for WIMPs interacting with heavy 
nuclei using data from the CDMS and DAMA experiments. The experimen¬ 
tal constraints (Higgs mass, CDM relic density) favour scenarios with large 
tan /?. For tan /? = 50 we found that the anapole moment of the lightest 
neutralino of the cMSSM has a value (P(10“^ — 10““^) fm which lies at 
the border of sensitivity of the current experimental searches and within the 
reach of the future experiments like XENONIT and LZ. 

The same kind of calculation can be performed in the context of other 
more complex models. This could be extremely valuable for discriminating 
not only between different interesting regions of the cMSSM, but also among 
different models. 

Thus, the anapole analysis could be useful to study the allowed parameter 
space of the cMSSM and other SUSY models. 
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Appendix: Scalar Three-point Function 


In this appendix, we analyse the Passarino-Veltman scalar three-point func¬ 
tion Co{q‘^ 


1^,93 which appears in the TDM calculation. 




^2 2 


Here denotes the photon transfered 4-momentum, x is the neutralino mass, 
and y and z are the masses of the particles running in the loop. 



Figure 9: Comparison between numerical (red line) and approximate (blue 
line) scalar three-point function Co(g^, ?/^), with x = 97.7 GeV, 

y = 415.4 GeV and z = 80.43 GeV. The analytical approximation (blue line) 
is only valid for —)■ 0. 

The corresponding plot for this Cq function can be seen in Fig. The 
red line shows the numerical solution, the blue line line represents the ap¬ 
proximate solution, i.e., the Taylor expansion around q^ = 0, which can be 
written as follows: 
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( 17 ) 


Co /) = Oo + + 0{q*). 

The coefficients a* are functions of the masses: 


tto — 


log (S) 


2 a ;2 


+ alogcc, 


(18) 


— y‘^x^ — 2z^x‘^ + 

Oil = - 

%x‘^z‘^{—x + y — z){x + y — z){—x + y + z){x + y + z) 

where 



12x4 


f 6 logo;, 

(19) 


+ iy^ — iz^ + ■\/—y^ + 2(x^ + z'^)y^ — (zP- — x'^)^ {ix^ — iy^ + iz'^ + \/—y"^ + 2{x’^ + z^)y’^ — (z^ — x^)^ 

^—ix^ + iy^ — iz'^ + ^—y'^ + 2(a;2 + z^)y‘^ — (z^ — x^)^ ^—ix^ — iy'^ + iz'^ + ■sZ—y'^ + 2(x'^ + z^)y^ — (z^ — x^)^ 

(20) 


i{x‘^ + y"^ — z^) 

2x2 a/— x4 + 2y'^x^ + 2z^x^ — y"^ — z"^ + 2y'^z^ 

and 


( 21 ) 


i(x^ + y’^ — z^){x*^ — Ay'^x'^ — 2z'^x'^ + — ‘iy'^z'^) 

12x'^{—x + y — z)(x + y — z)(—x + y + z)(x + y + z)^/—x'^ + 2y'^x‘^ + 2z‘^x'^ — y"^ — z'^ -f- 2y‘^z'^ 

( 22 ) 


References 

[1] A. M. Szelc, Acta Phys.Polon. B41, 1417 (2010), arXiv:1010.3918. 

[2] R. Schnee, (2011), arXiv:1101.5205. 

[3] K. Sigurdson, M. Doran, A. Kurylov, R. R. Caldwell, and 

M. Kamionkowski, Phys.Rev. D70, 083501 (2004), arXiv:astro- 

ph/0406355. 

[4] P. Ciafaloni et ai, JCAP 1103, 019 (2011), arXiv: 1009.0224. 


21 



[5] P. Ciafaloni and A. Urbano, Phys.Rev. D82, 043512 (2010), 
arXiv: 1001.3950. 

[6] P. Ciafaloni et ai, (2011), arXiv:1104.2996. 

[7] N. F. Bell, J. B. Dent, T. D. Jacques, and T. J. Weiler, Phys.Rev. D83, 
013001 (2011), arXiv:1009.2584. 

[8] N. F. Bell, J. B. Dent, T. D. Jacques, and T. J. Weiler, (2011), 
arXiv:1101.3357. 

[9] N. F. Bell et ai, (2011), arXiv: 1104.3823. 

[10] J. B. Dent, R. J. Scherrer, and T. J. Weiler, Phys.Rev. D78, 063509 
(2008), arXiv:0806.0370. 

[11] M. Kachelriess, P. Serpico, and M. Solberg, Phys.Rev. D80, 123533 
(2009), arXiv:0911.0001. 

[12] J. H. Heo, (2009), arXiv:0902.2643. 

[13] J. H. Heo, Phys.Lett. B693, 255 (2010), arXiv:0901.3815. 

[14] S. Liebler and W. Porod, Nucl.Phys. B849, 213 (2011), 

arXiv:1011.6163, * Temporary entry *. 

[15] V. Barger, W.-Y. Keung, and D. Marfatia, Phys.Lett. B696, 74 (2011), 
arXiv: 1007.4345. 

[16] V. Berezinsky, M. Kachelriess, and S. Ostapchenko, Phys.Rev.Lett. 89, 
171802 (2002), arXiv:hep-ph/0205218. 

[17] M. Kachelriess and P. Serpico, Phys.Rev. D76, 063516 (2007), 
arXiv:0707.0209. 

[18] L. Cabral-Rosetti, M. Mondragon, and E. Reyes Perez, AIP Conf.Proc. 
1116, 447 (2009). 

[19] L. Cabral-Rosetti, M. Mondragon, and E. Reyes Perez, J.Phys.Conf.Ser. 
315, 012007 (2011). 

[20] L. G. Cabral-Rosetti, M. Mondragon, and E. Reyes-Perez, (2012), 
arXiv: 1206.5052. 


22 



[21] L. Cabral-Rosetti, M. Mondragon, and E. Reyes Perez, J.Phys.Conf.Ser. 
485, 012019 (2014). 

[22] C. M. Ho and R. J. Scherrer, Phys.Lett. B722, 341 (2013), 
arXiv:1211.0503. 

[23] J. Kopp, L. Michaels, and J. Smirnov, JCAP 1404, 022 (2014), 
arXiv:1401.6457. 

[24] E. Del Nobile, G. B. Gelmini, P. Gondolo, and J.-H. Huh, (2014), 
arXiv: 1401.4508. 

[25] A. Ibarra, G. E. Yaguna, and O. Zapata, Phys. Rev. D93, 035012 
(2016), arXiv:1601.01163. 

[26] Y. Zeldovich, Zh. Eksp. Teor. Piz 33, 1531 (1957). 

[27] V. Dubovik and A. Gheshkov, Sov. J. Particles Nucl. 5, 318 (1975). 

[28] M. Abak and G. Aydin, Europhys. Lett. 4, 881 (1987). 

[29] V. M. Dubovik and V. E. Kuznetsov, Int. J. Mod. Phys. A13, 5257 
(1998), arXiv:hep-ph/9606258. 

[30] E. N. Bukina, V. M. Dubovik, and V. E. Kuznetsov, Phys.Lett. B435, 
134 (1998), arXiv:hep-ph/9805491. 

[31] B. K. Sahoo, T. Aoki, B. P. Das, and Y. Sakemi, (2015), 
arXiv: 1512.02055. 

[32] Y. V. Stadnik and V. V. Elambaum, Phys. Rev. D89, 043522 (2014), 
arXiv:1312.6667. 

[33] B. M. Roberts et al, Phys. Rev. D90, 096005 (2014), arXiv:1409.2564. 

[34] Y. V. Stadnik and V. V. Elambaum, Eur. Phys. J. C75, 110 (2015), 
arXiv: 1408.2184. 

[35] M. 1. Gresham and K. M. Zurek, Phys. Rev. D89, 123521 (2014), 
arXiv: 1401.3739. 

[36] E. O. Kamenetskii, Europhysics Letters 65, 269 (2004), arXiv:cond- 
mat:0305431. 


23 



[37] A. E. Miroshnichenko et ai, Nature communications 6 (2015). 

[38] C. H. Page, American Journal of Physics 39, 1039 (1971). 

[39] J. Corum, IEEE Ant. Prop. Int. Symp. 25, 832 (1987). 

[40] G. N. Afanasiev, Journal of Physics A Mathematical General 23, 5755 
(1990). 

[41] G. N. Afanasev and V. M. Dubovik, Sov. J. Part. Nucl. 23, 552 (1992). 

[42] N. J. Carron, American Journal of Physics 63, 717 (1995). 

[43] Y. Gao, C. M. Ho, and R. J. Scherrer, (2013), arXiv:1311.5630. 

[44] U. Haisch and F. Kahlhoefer, JCAP 1304, 050 (2013), arXiv:1302.4454. 

[45] D. Larson et ai, Astrophys.J.Suppl. 192, 16 (2011), arXiv: 1001.4635. 

[46] G. Jungman, M. Kamionkowski, and K. Griest, Phys.Rept. 267, 195 
(1996), arXiv:hep-ph/9506380. 

[47] K. A. Olive, p. 797 (2003), arXiv:astro-ph/0301505. 

[48] J. L. Feng, Ann.Rev.Astron.Astrophys. 48, 495 (2010), arXiv:1003.0904. 

[49] S. P. Martin, (1997), arXiv:hep-ph/9709356. 

[50] I. J. Aitchison, (2005), arXiv:hep-ph/0505105, Notes of Lectures for 
Graduate Students in Particle Physics Oxford, 2004 and 2005. 

[51] L. Roszkowski, E. M. Sessolo, and A. J. Williams, JHEP 08, 067 (2014), 
arXiv: 1405.4289. 

[52] A. Fowlie, Eur. Phys. J. C74, 3105 (2014), arXiv: 1407.7534. 

[53] C. Balazs, A. Buckley, D. Carter, B. Farmer, and M. White, Eur. Phys. 
J. C73, 2563 (2013), arXiv: 1205.1568. 

[54] C. Strege et ai, JCAP 1304, 013 (2013), arXiv:1212.2636. 

[55] M. E. Cabrera, J. A. Casas, and R. R. de Austri, JHEP 07, 182 (2013), 
arXiv:1212.4821. 


24 



[56] O. Buchmueller et al, Eur. Phys. J. C72, 2243 (2012), arXiv:1207.7315. 

[57] P. Bechtle et al, How alive is constrained SUSY really?, in 37th Inter¬ 
national Conference on High Energy Physics (ICHEP 2014) Valencia, 
Spain, July 2-9, 2014, 2014, arXiv: 1410.6035. 

[58] E. A. Bagnaschi et al, Eur. Phys. J. C75, 500 (2015), arXiv:1508.01173. 

[59] P. Bechtle et al, JHEP 06, 098 (2012), arXiv:1204.4199. 

[60] C. Strege et al, JCAP 1203, 030 (2012), arXiv:1112.4192. 

[61] W. Abdallah and S. Khalil, (2015), arXiv: 1509.07031. 

[62] M. E. Cabrera-Catalan, S. Ando, C. Weniger, and E. Zandanel, Phys. 
Rev. D92, 035018 (2015), arXiv:1503.00599. 

[63] J. Ellis et al, (2015), arXiv:1509.08838. 

[64] B. C. Allanach, D. P. George, and B. Gripaios, JHEP 07, 098 (2013), 
arXiv: 1304.5462. 

[65] B. G. Allanach, D. P. George, and B. Nachman, JHEP 02, 031 (2014), 
arXiv:1311.3960. 

[66] O. Buchmueller et al., Eur. Phys. J. C74, 2809 (2014), arXiv:1312.5233. 

[67] J. Ellis, Prospects for Supersymmetry at the LHG & Beyond, in 
18th International Conference Prom the Planck Scale to the Elec- 
troweak Scale (Planck 2015) loannina, Creece, May 25-29, 2015, 2015, 
arXiv: 1510.06204. 

[68] E. N. Bukina, V. M. Dubovik, and V. E. Kuznetsov, JINR-E4-98-95. 

[69] A. Gongora and E. Ley-Koo, Rev. Mex. Eis. E 52, 188 (2006). 

[70] M. Pospelov and T. ter Veldhuis, Phys. Lett. B480, 181 (2000), 
arXiv:hep-ph/0003010. 

[71] R. Bernabei et al, Phys. Lett. B389, 757 (1996). 

[72] GDMS, R. Abusaidi et al, Phys. Rev. Lett. 84, 5699 (2000), arXiv:astro- 
ph/0002471. 


25 



[73] M. Klasen, M. Pohl, and G. Sigl, Prog. Part. Nucl. Phys. 85, 1 (2015), 
arXiv: 1507.03800. 

[74] XENONIOO, E. Aprile et al, Phys. Rev. Lett. 109, 181301 (2012), 
arXiv: 1207.5988. 

[75] LUX, D. S. Akerib et al, Phys. Rev. Lett. 112, 091303 (2014), 
arXiv:1310.8214. 

[76] XENON, F. A. Davide, EPJ Web Conf. 95, 04019 (2015). 

[77] LZ, D. S. Akerib et al, (2015), arXiv: 1509.02910. 

[78] XENONIOO, G. Kessler, XENONIOO and XENONIT Dark Matter 
Search with Liquid Xenon, in Proceedings, 20th International Confer¬ 
ence on Particles and Nuclei (PANIC If), pp. 357-360, 2014. 

[79] LUX, J. E. Y. Dobson, Searching for Dark Matter with the LUX exper¬ 
iment, in Proceedings, 20th International Conference on Particles and 
Nuclei (PANIC If), pp. 373-377, 2014. 

[80] R. Mertig, M. Rohm, and A. Denner, Gomput.Phys.Gommun. 64, 345 
(1991). 

[81] T. Hahn and M. Perez-Victoria, Gomput. Phys. Gommun. 118, 153 
(1999), arXiv:hep-ph/9807565. 

[82] A. Djouadi, J.-L. Kneur, and G. Moultaka, Gomput. Phys. Gommun. 
176, 426 (2007), arXiv:hep-ph/0211331. 

[83] ATLAS Gollaboration, G. Aad et al, Phys.Lett. B716, 1 (2012), 
arXiv: 1207.7214. 

[84] ATLAS Gollaboration, (2013), ATLAS-GONF-2013-014, ATLAS-GOM- 
GONF-2013-025. 

[85] GMS Gollaboration, S. Ghatrchyan et al, Phys.Lett. B716, 30 (2012), 
arXiv: 1207.7235. 

[86] GMS Gollaboration, S. Ghatrchyan et al, (2013), arXiv:1303.4571. 

[87] GMS and LHGb Gollaborations, (2013), 

http://cdsweb.cern.ch/record/1374913/files/BPH-11-019-pas.pdf. 


26 



[88] A. Arbey, M. Battaglia, F. Mahmoudi, and D. Martinez Santos, 
Phys.Rev. D87, 035026 (2013), arXiv: 1212.4887. 

[89] (2012), Heavy Flavour Averaging Group, see: 

WWW. slac.stanford.edu/xorg/hfag/. 

[90] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, Gomput. Phys. 
Gommun. 176, 367 (2007), arXiv:hep-ph/0607059. 

[91] A. Fowlie et ai, Phys. Rev. D86, 075010 (2012), arXiv: 1206.0264. 

[92] L. G. Gabral-Rosetti and M. A. Sanchiz-Lozano, J.Gomput.Appl.Math. 
115, 93 (2000), arXiv:hep-ph/9809213. 

[93] L. G. Gabral-Rosetti and M. A. Sanchiz-Lozano, J.Phys.Gonf.Ser. 37, 
82 (2006), arXiv:hep-ph/0206081. 


27 



